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The resistive transition of granular high-Tc superconductors, characterized by either weak 
(YBCO-like) or strong (MgB2-hke) hnks, occurs through a series of avalanche- type current den- 
sity rearrangements. These rearrangements correspond to the creation of resistive layers, crossing 
the whole specimen approximately orthogonal to the current density direction, due to the simul- 
taneous transition of a large number of weak-links or grains. The present work shows that exact 
solution of the Kirchhoff equations for strongly and weakly linked networks of nonlinear resistors, 
with Josephson junction characteristics, yield the subsequent formation of resistive layers within 
the superconductive matrix as temperature increases. Furthermore, the voltage noise observed at 
the transition is related to the resistive layer formation process. The noise intensity is estimated 
from the superposition of voltage drop elementary events related to the subsequent resistive layers. 
At the end of the transition, the layers mix-up, the step amplitude decreases and the resistance 
curve smoothes. This results in the suppression of noise, as experimentally found. Remarkably, a 
scaling law for the noise intensity with the network size is argued. It allows to extend the results to 
networks with arbitrary size and, thus, to real specimens. 

PACS numbers: 74.40.-|-k, 74.78.Bz, 74.81.-g 



I. INTRODUCTION 

The supcrconductivc-normal state transition may oc- 
cur according to diverse mechanisms, depending on phys- 
ical conditions, material type and structure. In type II 
superconductors at T <C Tc, where is the critical tem- 
perature in the absence of magnetic field and current, 
the transition occurs when fluxoids, injected by external 
magnetic fields or strong bias current densities, begin to 
move causing energy losses and heating. This is rele- 
vant for the development of high-field superconducting 
magnetsii^i^ii^i^i^i^ii^ii^. A different transition mechanism 
may occur when temperature is close to at low cur- 
rent density. In this case, an intermediate state may be 
obtained, characterized by a mixture of superconductive 
and normal domains. This situation was first studied by 
Landau and Ginzburg in metalsi^iii. Recently, it has 
been considered to explain the excess noise in metallic or 
high-Tc superconductors transition edge sensors (TES) 
used as bolometers to detect electromagnetic radiation 
at the level of single photon a^^'^'^i^^ . 
The excess noise observed during a transition sheds 
light on the microscopic processes underlying the tran- 
sition itselfi^ii^iiiii^. Ini^, the noise observed during the 
superconductor-normal transition in granular MgB2 films 
has been ascribed to the subsequent formation of resistive 
layers, with grains in the normal or in the intermediate 
state, between equipotential superconducting domains. 
The excess noise derives from the fact that each elemen- 
tary event -the formation of a layer- implies the simul- 
taneous resistive transition of several grains and, thus, 
gives rise to a voltage pulse of rather high amplitude 



[avalanche noise). 

The present work is addressed to simulate the transi- 
tion events occurring at granular level responsible for 
the avalanche-type noise in YBCO-like and MgB2-like 
supcrconductorsi^iiS. The superconducting material is 
modeled as a network of nonlinear resistors, having 
Josephson junction current-voltage (/ — V) characteris- 
tics with gaussian distribution of critical currents. The 
nonlinear resistors represent either weak links between 
grains (YBCO-like) or grains with strong links (MgB2- 
like). In the strong link case, a couple or triple of re- 
sistors is used to represent two or three current compo- 
nents flowing through each grain respectively for two- 
(2D) and three-dimensional (3D) networks. The solu- 
tions of Kirchhoff equations for these networks are found 
by an iterative routine described in the next section. The 
main results of this analysis are: 

1. the resistive transition undergoes discrete step-like 
increments both in weak and strong link materials. 
The steps correspond to the creation of resistive 
layers constituted by grains or weak links in the 
normal or in the intermediate state. As temper- 
ature increases, grains or weak-links in the inter- 
mediate state gradually switch to the normal state. 
The trailing edge of the resistive transition grows 
more smoothly in MgB2-like than in YBCO-like 
networks. This fact is related to the higher cor- 
relation when the elementary transition events oc- 
cur in triplets rather than in independent nonlinear 
resistors. 

2. the abrupt formation of resistive layers causes the 
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large voltage noise observed at the transition in 
these materials. At the end of the transition, the re- 
sistive layers mix up. The resistance steps become 
smaller and the transition curve smoother. This 
smoothing results in noise suppression. A scaling 
law for the noise intensity is proposed in order to 
extend the results to larger networks, representing 
real materials. This effect was simply assumed ini^. 
Here it is shown that the transition noise can be 
estimated once grain size and critical current dis- 
tribution are defined. 



II. NETWORKS OF STRONG AND WEAK 
LINKS 

Before describing the details of the simulations, we 
provide a description of the main physical parame- 
ters relevant to the electronic behavior of granular 
superconductors22i2i. In particular, it is worthy to re- 
mind that the phase transition boundary of granular su- 
perconductors is set by the value of the dimensionlcss 
tunneling conductance g: 
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where G is the average tunneling conductance between 
adjacent grains and /h the quantum conductance. 
Experiments show that samples with the normal state 
conductance greater than the quantum conductance 
(i.e. with the g ^ 1) become superconducting at low 
temperature^^, regardless of the ratio of Josephson J and 
Coulomb Ec energies, defined respectively as: 



with A the superconductor gap, and 



Ec = — 
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with Cj the grain capacitance. This phenomenon can 
be accounted by the electron tunneling between grains, 
in addition to the Josephson coupling2^. The additional 
dissipative tunneling channel results in a reduction of the 
Coulomb energy to: 
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known as the effective Coulomb energy of the grain. By 
comparing the Eqs. ^ and one can notice that for 
5^1, J is always larger than Ec, implying a supercon- 
ducting ground state, regardless of the Coulomb energy 
Ec- For g ^ 1, the granular superconductor can then 
be modeled within the mean- field BCS theory. Thus, its 
critical temperature is approximately given by the sin- 
gle grain BCS critical temperature 2A = 3.53fcTc. Con- 
versely, for g ^ 1, the phase transition boundary between 



insulating and superconducting states is controlled by the 
ratio between J and Ec- In this condition, by using a 
mean field approach, the critical temperature is given by 
Tc = (l/4)z7r(7A, with z the coordination number of the 
lattice^iSi. 

The superconductor-normal transition in thin granular 
films with g ^ 1 can be modeled in terms of rcsistivcly 
shunted Josephson junctions, whose state is controlled 
only by the value of the normal resistance, rather than 
by the Coulomb and Josephson energies. The simulations 
presented in this work have been performed in the regime 
g ^ 1, to guarantee the onset of a superconductivity 
state at low temperature. 
In order to simulate the superconductor-normal transi- 
tion in granular materials, two types of networks, shown 
in Figs.[lja) and (b), are considered. The networks are 
constituted by non-linear resistors, with Josephson junc- 
tion characteristics, biased by a constant current gener- 
ator. The resistive transition is estimated by solving a 
system of Kirchhoff equations, at varying temperature, 
for each network. 

The network of Fig. [Ija) refers to YBCO-like materials, 
characterized by weak-links^^. In these materials, the 
transition occurs in two separated steps: first, at lower 
temperatures, for the weak-links and, then, at higher 
temperatures, for the grains. The network of Fig. [BJa) 
is used to model the first stage of the transition, which 
involves only the weak-links, while the grains remain su- 
perconductive. 

The network of Fig. mb) refers to MgB2-like 
materials2^i^i21, whose transition involves directly the 
grains. Each triplet of resistors, outlined by the red cir- 
cle, represents a grain. Since the current density within 
the grain may have any direction, the three resistors give 
a basis of three components of the current density for 
each grain. The current density sets the state (super- 
conductive, intermediate, normal) of the grain according 
to its / — 1^ characteristic. For the sake of simplicity, 
all the grains arc assumed to be isotropic and with the 
same average size, therefore the anisotropy effects are 
disregarded22i22ii2£. This assumption is not limitative for 
what concerns the main aspects of the transition. It al- 
lows one to define a critical current Il{T) characterizing 
the grain i, according to a Gaussian distribution, and a 
normal state resistance Rg equal for all the grains. In 
real specimens, small changes of the grain stoichiome- 
try influence the critical current more than the normal 
state resistivity. The spread of the distribution of the 
critical currents and temperatures is responsible of the 
slope of transition curve^i. The normal state resistance 
Ro is achieved when the current crossing the grain or 
the weak- link exceeds Il{T). The intermediate states are 
characterized by current Il.{T) and voltage drop between 
and V^iT). The I — V characteristic of each non-linear 
resistor, representing a grain or a weak-link, is completely 
defined by the quantities: /*(T) and R. The quantity 
Vc = Il{T)R is directly related to the Josephson time 
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(b) 



FIG. 1: Scheme of 3D networks of granular superconductors 
with: (a) weak-hnks (YBCO-like) and (b) strong-hnks, where 
the transition involves directly the grains (MgB2-like) . The 
networks contain 48 grains. In (a), the grains are assumed 
to remain in the superconducting state during the transition, 
and correspond to the nodes of the network. Each link be- 
tween grains is a nonlinear resistor with the I—V characteris- 
tic represented in Fig[2]In (b), the I — V characteristic given 
in Fig[2] concerns the whole grain, which is represented by 
a triplet of orthogonal resistors (outlined by the red circle). 
The electrical conductance of these resistors is determined 
from the I — V characteristic of the grain by calculating the 
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potential drop across the grain i asVi = | 



where 



j identifies the three grains linked to the grain i. Since the 
grain is assumed to be isotropic, the conductances of the 3 
resistors are assumed equal. The first (1) and last node (50) 
correspond to the electrodes. 



constant by: 
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for the intermediate states (0 < u < Vc) and 
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for the normal states {v > Vc). These time constants 
define the characteristic switching time during the tran- 
sition, and thus, are ultimately related to the behavior 
of noise. In the next section, the resistive transition has 




FIG. 2: I — V characteristics used for the simulation of the 
transition in the grain o weak-link networks. The resistance 
is assumed to be negligible (lO'^'^il), when I < Il{T), while 
for I > IliT) the normal state resistance is assumed to be 
Ro, equal for all the grains or links. The intermediate states 
occur for voltage drops between and V^{T) — RoIl{T) and 
current I — Il{T). The critical current IliT) is distributed 
according to a Gaussian distribution function for the difi^er- 
ent grains. The transition from the superconductive to the 
intermediate state is set by the value of 7c. The transition 
from the intermediate to the normal states is determined by 
the product Vc = RIc- Therefore the Gaussian distribution 
of the Ic is enough to ensure the randomness of the product 
Vc = RIc for all the grains. Furthermore, the quantity Vc 
is directly related to the Josephson time constant tj, given 
by the Eqs. ([Sjl and (|6]), that define the elementary switching 
time of the transition, and thus, are ultimately related to the 
behavior of noise. 



been simulated in networks with (a) underdamped, (b) 
overdamped and (c) general I — V characteristics, that 
are characterized by the Stewart-McCumber parameter 
13c = tuc/tj, where t^q a-nd tj are the capacitance and 
Josephson time constant. /3c ^ 1 in the case (a), /3c ^ 1 
in the case (b) and /?c 1 in the case (c). In particular, 
the onset of hysteresis has been analyzed upon cooling 
the granular system from the normal to the supercon- 
ductive state. 

When the transition involves the grains (strong- links), 
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where 



the current is given by: li 

corresponds to the current flowing from the grain i to its 
neighboring grains j through each resistor of the triplet 
(see Fig. [T](b)). The I — V characteristic is then used to 
find the value of the three resistors by means of an itera- 
tive routine to solve KirchhofF equations. The grains are 
assumed isotropic, thus the three resistors representing 
the grain will have identical I — V characteristics. 
The simulations are carried on at constant bias current. 
The transition is caused by the temperature increase, 
which reduces the critical currents of the grains or weak- 
links according to the following linearized equation: 



IciT) = Pc, 



T 

1 

T, 



(7) 



where I^^ is the low-temperature critical current, dis- 
tributed according to a Gaussian function with standard 
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The preliminary steps of the simulations are as follows: 



1 . the list of all the No nodes of the network is created. 

2. the Gaussian distribution for the critical current 
J^o is introduced. In Matlab, the vector I^^ is de- 
fined by: /*Q = /co (1 + randn{N, 1)), where N is 
the number of junctions (individual resistors) for 
the network (a) or the number of grains (triplets 
of resistors) for the network (b). The quantity 
randn{N, 1) defines a set of N random numbers ex- 
tracted from a gaussian distribution having mean 
value and variance 1. 

Then the iterative calculations are implemented as fol- 
lows: 

• the vector Wo of the tentative potential values is 
defined for all the No nodes. 



for the network of Fig[lja), by using the I — V 
characteristics, a conductance value Gy for each 
resistor between the nodes i and j is calculated 



else, for the network of Fig[l] (b), the conductance 
values, common to the three resistors representing 
each grain i, are calculated from the I ~V charac- 
teristics, by using the voltage drop: 
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Once the Gij are known, the entries of the conductance 
matrix G are: 



Gi: 



—Gij contiguous) , 

not contiguous) 

Eg. 



(9a) 
(9b) 

(9c) 



where Gik are the conductances of the resistors connected 
to the node i. 

Then, a new vector of node potentials Wi is evaluated 
by solving the equation: 



the range 10~^ < Smin < 10""'^"'^ to check that the value of 
Emin did not appreciably change the final solution. The 
total network resistance R is then given by Wn{l)/Ib for 
each value of T/Tc, where W„(l) is the potential drop at 
the contact ends. 

The potential drops at the ends of each resistor for the 
case (a) and across the grain for the case (b) are com- 
pared to the values of the potential in the corresponding 
I — V characteristics. Therefore, it is possible to distin- 
guish weak-links or grains respectively in the supercon- 
ducting, normal or intermediate state. Before discussing 
the simulations results, it is worthy to point to the dif- 
ferent behavior of the two networks by introducing the 
intragrain conductance gintr- For standard granular sys- 
tem, the condition g <C gintr holds. 

The intragrain conductance of the weak-link network 
shown in Fig. 1(a) is much greater than 1 {gintr ^ !)• 
The intragrain region is indeed assumed to remain in the 
superconducting state, since the transition occurs only at 
the weak-links. Conversely for the strong-link network of 
Fig. 1(b) the condition g ~ gintr holds, corresponding 
to an homogenously disordered granular system. This 
condition is consistent with the electronic properties of 
MgB2-like superconductors^^. The intragrain conduc- 
tance gintr is related to the single grain Thouless energy 
Eth and to the interlevel spacing S through 



gintr = Exh/ S 



(11) 



G-Wi= /„ 



(10) 



When the energy Exh exceeds the mean level spacing 
6, it is gintr ^ 1- The Thouless energy is defined by 
Eth — Do/a? , with Do and a the diffusion coefhcient and 
the radius of the grain. The interlevel spacing is defined 
5 = l/{vV) with V and V the density of states at the 
Fermi energy and the volume of the grain. The intragrain 
conductance strongly depends on the dirtiness of mate- 
rial and the radius of grain. These aspects are indeed 
relevant for MgB2-like materials whose critical temper- 
ature is strongly dependent on material quality, atomic 
radii and cell siz o^^i^^ . 



III. RESULTS 



with respect to Wi. linj is a vector of dimension No, 
whose elements are zero except the first one equal to the 
bias current /{,. It represents the external current in- 
jected into node 1. The last node is grounded. 
Then, the new set of potentials Wi allows to evalu- 
ate a new set of Gij and a new conductance matrix 
G. From Eq. ([TU)l an updated vector 14^2 is then ob- 
tained. The iteration is repeated until the quantity 
e = \Wn ~ Wn-i\/\Wn\ bccomes smaller than a value 
Emin chosen to exit from the loop. In the present work, 
the simulations have been performed by varying £mm in 



Here, the successive stages of the resistive transition 
are simulated in granular superconducting materials ei- 
ther with strong or with weak links. The supercon- 
ducting material is represented as a network of nonlin- 
ear resistors having resistively and capacitively shunted 
Josephson junction characteristic9^2iSi2i. We report the 
results of different simulations, carried on with 2D and 
3D networks, both for grain and weak-link transition. In 
the simulations, the transition occurs by increasing the 
temperature, in proximity of the critical temperature T^,, 
starting from the superconductive state. 
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A. Resistive layers in strong and weak-links 
networks 

Figs. [3] and m refer to the resistive transition of a two- 
dimensional 30 X 30 network, representing a granular su- 
perconducting film of 900 grains characterized by strong 
links (MgB2 type). 

Figs. [5] and [5] refer to a two-dimensional 30 x 30 network, 
representing a superconducting film of 900 grains char- 
acterized by weak-links (YBCO type). 
R and T are expressed as reduced quantities, namely as 
R/Ro, T/Tc- The relevant energy values and the param- 
eters used for the simulations are reported in Table [H in 
Table HIl and / or in the figure captions. 

In Fig. m at the beginning of the transition, the net- 
work resistance is zero, since all the grains are in the 
superconductive state. By effect of the temperature in- 
crease, a layer of grains either in the normal resistive 
(blue) or in the intermediate (green) state, crossing the 
whole film, is generated (Fig. [3] (a)). This layer must 
separate two equipotential superconductive regions and, 
thus, the potential drop must be constant along the layer. 
Since the grains have different critical currents, the layer 
starts to form when the sum of critical currents of its 
grains equals the bias current. The grain (or the weak- 
link) with the lowest critical current becomes resistive 
and set the voltage drop of the other grains in the layer. 
As temperature increases and critical current decreases, 
more and more grains in the intermediate state gradu- 
ally switch to the resistive states and the layer resistance 
increases. 

As shown in Fig. [3l a resistive layer contains at least 
one resistive (blue) dot and many intermediate (green) 
dots. Superconductive (red) dots are excluded since they 
would constitute a short. The formation of a resistive 
layer corresponds to a step in the R vs T curve, as it 
can be seen by comparing Figs. [3] and SI Upon further 
temperature increase, other layers are created until the 
whole film undergoes the transition to the normal state. 
At the beginning of the transition the layers are well 
separated and have a thickness of approximately one 
grain. Correspondingly, the resistance steps shown in 
Fig. |4] obey, as a good approximation, to a scaling law 
{R/Ro = 1/30 in the present case). At the transition 
end, there is an intricacy of different layers and the re- 
sistance increases smoothly with the temperature. 
Figs. [5] and [6] correspond to granular superconductors 
characterized by weak-links (YBCO-like). The simula- 
tion refers to the resistive transition of the weak-links. 
The grains, represented by the nodes of the network, re- 
main in the superconductive state. Also in this case the 
transition occurs through the formation of resistive lay- 
ers corresponding to resistance steps in the Rvs T curve. 
Figs. [7] and [5] report simulations carried onlOxIOxlO 
3D networks, representing superconductor films of 1000 
grains respectively with strong and weak links. The pres- 
ence of about 10 steps is expected from the scaling law 
holding before mixing up of the layers {R/Ro = 1/10). 
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FIG. 3: Different stages of the superconductor-normal transi- 
tion in a 2-dimensional MgB2-like network with 30 x 30 grains. 
Red dots represent the superconductive grains, blue dots rep- 
resent the resistive grains, green dots represent grains in the 
intermediate state. In (a), the first resistive layer (a strip in 
2D) is formed, which corresponds to the first step in the R vs 
T curve of Fig. |4l In (b) , at slightly higher temperature the 
appearance of the second layer is shown, corresponding to the 
second step in Fig.|3](a) etc. (c) The formation of more layers 
is shown. In (d) , the situation at the transition end is shown. 
The layers mix-up and the resistance steps become smoother, 
as shown in Fig.|4](b). The parameters used in the simulation 
are : Tc = Z'dK; Ro = 0.32^; = lmA\ ho = 1.7mA. 
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TABLE I: Relevant energy scales. The value of the capacitance used for the calculation of the Coulomb energy are 2.7-10 ''pF < 
Cj < 2.7 ■ IQ-^pF. 





Y BOO 


MgB2 






TT bands 


(7 bands 


Superconductive gap (A) 


10 meV ^ 20 meV 


1.2meV 3.7 meV 


6.4 ~ 7.2 meV 


Critical temperature (Tc) 


65.8 K ^ 131.5 K 


7.9 K^24.3K 


42.1 K^47.4 K 


Ambegaokar-Baratoff product (14 = 7rA/(2e)) 


15.7 mV ^ 31.4 mV 


1.9 mV^ 5.8 mV 


10.1 mV ^ 11.3 mV 


Coulomb energy (_Ec) 


0.1 eV ^ 0.1 • 10"* eV 


0.1 eV ~ 0.1 ■ 10"* eV 


0.1 eV ^ 0.1 ■ 10"'' eV 


Effective Coulomb energy (Ec) 


3.0 /ieV ^ 5.0 ^eV 


0.03 fieV^ 0.09 ^eV 


0.16 ^eV ^ 0.18 neV 


Josephson coupling energy (J) 


63.3 eV ~ 126.6 eV 


76.0 eV-^ 234.2 eV 


405.2 eV ^ 455.8 eV 



TABLE 11: Simulation parameters. 





min 


YBCO 

max 


MgB2 

min 


max 


Low temperature Critical current (/co) 


1.0 mA 


10 mA 


1.0 mA 


10 mA 


Normal state resistance (Ro) 


0.1 Q, 


3.0 n 


0.1 n 


i.on 


Critical voltage (K = Roico) 


0.1 mV 


10 mV 


0.1 mV 


10 mV 


Dimensionless tunneling conductance (g) 


1.3 ■ 10^ 


4.3 • 10^ 


1.3 ■ 10^ 


1.3 ■ 10" 


Critical temperature (Tc) 


65.8 K 


131.5 K 


11.8 K 


44.7 K 




FIG. 4: Resistive transition of the two-dimensional network 
of 30 X 30 superconducting grains shown in Fig. |3] Figs. |4] 
(a) and (b) represent a zoom of the curve in two temperature 
intervals at the beginning and near the end of the transition. 
It may be noticed that the values of T/Tc in correspondence of 
the steps in Fig.|3](a) correspond also to the layers represented 
in Figs. [3](a),|3](b),|3](c). 



B. Hysteresis effects 

So far, we have been concerned with the 
superconductive-resistive transition as temperature 
increases with the main aim to investigate the successive 
formation of layers. Here, wc investigate what the 
algorithm can predict when temperature is lowered and 
the superconductive final state is achieved starting from 



the normal one, thus addressing the hysteresis onset. 
For this purpose, it is necessary to distinguish the I ~ V 
characteristics of rcsistivcly shunted undcrdamped (a), 
overdamped (b) and generalized (c) Josephson junctions 
(shown respectively Figs. Eta), (b),(c))^>23i34. The curve 
I — V (a) is hysteretic, the curve (b) shows no hysteresis, 
while the curve (c) exhibits partial hysteresis. 
We have routinely solved the Kirchhoff equations 
of the strong and weak links networks by using the 
underdamped, overdamped and generalized I ~ V 
characteristics and implementing a heating-cooling cycle 
around the critical temperature Tc- For all the three 
cases: (i) the conductance is G = lO^^S at I < Il{T), 
{ii) the normal state conductance Go = l/i?o at 
I > Il{T) has been varied in the range reported in table 
im (Hi) G and Go are much greater than the quantum 
conductance (i.e. .g 3> 1 always). 

For the underdamped I ~ V characteristics (a), the 
intermediate states are characterized by voltage drop 
in the range < V < Vl{T) and current equal to 
/ = Ic{T). The intermediate states correspond to the 
coexistence of superconducting and normal domains. 
Upon current (voltage) decrease starting from the 
normal state, the behavior is always normal resistive, 
implying that the system reaches the superconductive 
ground state without exploring intermediate states. 
For the overdamped I — V characteristics (b), the 
intermediate states are characterized by voltage drop 
in the range < V < 2V^{T) and current in the range 
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FIG. 6: Same as Fig. |4] but for the network of weak-links 
shown in Fig. (5] 



2 



(b) /f 



(a) 



(b) 



(a) 



0.30 
0.20 
0.10 



0.9800 9801 9802 0.9803 



10 



15 



20 



25 



30 



0.00 - 

0.9796 0.9800 0.9804 0.9808 

T/Tc 



FIG. 7: Resistive transition of a three-dimensional network of 
10 X 10 X 10 grains. The parameters used in this simulation 
are Ro ~ 0.32f2 and Ico ~ 1.7mA. Again each step should 
correspond to the creation of a layer of grains either in the 
normal or in the intermediate state through the network cross- 
section. 



FIG. 5: Different stages of the superconductor-normal tran- 
sition in a 2-dimensional YBCO-like network with 30 x 30 
grains. Red lines represent weak-links in the superconduct- 
ing state, blue lines represent resistive weak-links, while green 
lines represent weak-links in the intermediate state, (a) shows 
the formation of the first resistive layer (a strip in 2D) and 
corresponds to the first step in the RvsT curve of Fig. |6] (a) . 
(b) and (c) show the transition at a temperature immedi- 
ately following the appearance of the second and third layer 
and corresponds to the second and third step in Fig. |6] (a) 
etc. (d) shows the transition end, where the layers become 
mixed up and the resistance steps become smaller, as shown 
in Fig. |6] (b) . The parameters used in the calculations are: 
Tc = 65K; Ro = 0.32^; h = 1mA; ho = 2.1mA. 



O 0.60 

2 



(b) ; 



(a) 



(b) 



0.70 
0.60 



9847 0.9847 0.9i 



(a) 



3.9844 0.9846 0.9848 0.9850 0.9852 

T/Tc 



. 0.9844 0.9845 0.9846 



FIG. 8: Same as Fig. [7] but for weak-links with the parame- 
ters: Ro = 0.32f2 and Ico ~ 1.7mA. 
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li{T) < I < li{2V^{T)), as described by the function: 



V = IR^l-(^^J , (12) 

instead of a constant value. The behavior of the over- 
damped Josephson junction is the same upon increasing 
and decreasing the current (voltage). 
Fig. M^c) corresponds to the general case, the I — V 
curve is partly hysteretic. Upon heating, the intermedi- 
ate states are characterized by a voltage drop in the range 
< V" < 2V^{T) and current equal to /^(T). Conversely, 
upon cooling the intermediate states are described by Eq. 

m- 

Fig. [To] shows the resistive transition during a heating- 
cooling cycle in the case of a 2-dimcnsional network with 
weak links. In particular. Fig. (TU] (a) refers to a network 
with underdamped weak links, where the maximum hys- 
teresis effect can be observed. Fig. [10] (b) refers to a 
network with overdamped weak links, and no hysteresis 
is observed. Fig. [TU] (c) refers to a network of weak links 
with generalized Josephson junction characteristic, where 
the amount of hysteresis is an average of the previous two 
cases. 



C. Avalanche Noise at the Transition 

The described results show that on increasing the tem- 
perature from the superconductive state, subsequent re- 
sistive layers are formed, until the whole specimen be- 
comes normal. These layers are abruptly formed across 
the networks and correspond to step-like resistance in- 
crements in the i? vs T curves. Each resistance step 
involves the simultaneous transition of all the grains in 
a layer, and corresponds to a voltage pulse at the end of 
the network, when a constant bias current lb is applied. 
The number of these pulses during the transition is of 
the order of the number of grains or weak-links along the 
current direction. The large transition noise is due to the 
random superposition of these voltage pulses. 
For a given value of the normal state resistance, the step 
resistance and the step voltage amplitude are inversely 
proportional to the number of steps. In real specimens 
with huge number of grains, the resistance steps can- 
not be resolved by static measurements of the transition 
curve. Conversely, the noise is a measure of the tran- 
sition dynamics at granular level. The noise amplitude 
depends on the number and amplitude of discrete volt- 
age steps and thus permits to justify the step-transition 
model. Assuming a Poisson distribution of the pulses, the 
power spectrum of the noise is given by (Campbell 
theorem) : 



^{u;) = i^\S^{iu)\^ , (13) 

where v represents the average number of pulses per unit 
time, and |S't,(a;)|2 the average square modulus of the 




|b) 



FIG. 9: Josephson junction I — V characteristics for grains 
or weak-links in case of underdamped (a), overdamped (b), 
generalized (c). I^in depends on the Stewart- McCumber pa- 
rameter Pc and ranges from II and for /3c > 0, where 
/3c = trc/tj, where trc and tj are the capacitance and 
Josephson time constant respectively. /3c ^ 1 in the case (a), 
/3c ^ 1 in the case (b) and /?c ~ 1 in the case (c) 



Fourier transform of each pulse. 

Real materials correspond to very large networks, whose 
number of nodes is obtained by dividing the specimen 
dimensions by the average grain size. The amplitude of 
the voltage pulses is inversely proportional to the number 
of steps, while ly is directly proportional to the number 
of steps along the transition curve. Since, in the low 
frequency limit (w ^ 0), [^^(O)!^ is proportional to the 
square amplitude of the voltage pulse, it turns out that 
the noise amplitude is inversely proportional to the num- 
ber of steps. This shows that, for a given value of the 
network resistance, the voltage noise amplitude in the 
limit of low-frequency is inversely proportional the num- 
ber of grains along the direction of the flowing current. 
This, in other words, means that superconductors with 
smaller grains are characterized by a lower intensity of 
the transition noise. 
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FIG. 10: Resistive transition of the two-dimensional network 
of 30 X 30 superconducting weak-links when temperature in- 
creases and decreases. The degree of disorder is 0.1 like the 
simulations. The external current I is kept constant at 1mA. 
The Kirchhoff equation iteration is implemented with (a) un- 
derdamped Josephson junctions [Pc S> 1), (b) overdamped 
Josephson junctions [Pc <^ 1) and general Josephson junc- 
tions {Pc ~ 1). One can note respectively hysteresis (a), no 
hysteresis (b) and partial hysteresis upon cooling the material 
from the normal to the superconductive state. 



tween the layer formation and the transition noise at 
low- frequencies. The present approach gives exact nu- 
merical solutions for the transition. In addition it allows 
to evidence the decrease of noise towards the transition 
end. By representing the superconducting film as a net- 
work of non linear resistors, it is possible to evaluate 
how the noisiness decreases at the end of the resistive 
transition, according to the variance of the distribution 
of the grain or of the weak-links critical currents. This 
is a crucial issue for the development of superconductor 
based sensor a'^^i'^^'^^ . The steepness of the R vs T curve 
gives higher photon detection signals (photoresponse) at 
the expenses of an increase of the noise. Moreover, the 
resistance steps, corresponding to each layer formation, 
are visibly more squared and sharp for weak links than 
for strong-links transition. This fact may be related to 
the slope of the relative voltage noise spectra reported in 
RefsjiSiii^. The power spectra are and l//^-sloped 
in the range between a few Hz and IkHz respectively for 
MgB2 and YBCO. Since the power spectrum of a random 
staircase signal (i.e. a sequence of Poisson distributed ex- 
ponential pulses, whose time decay tends to oo) is exactly 
l//^-sloped (i.e. a lorentzian function whose cut-off fre- 
quency tends to 0), the rounding of the pulse trailing 
edge produces a steeper decay of the power spectrum, 
which tends to the l//'^-slope. As a conclusion, it may 
be stated that the representation of granular supercon- 
ductors as a network of nonlinear resistors with resistively 
shunted Josephson junction characteristics add clues to 
the dynamics of the transition process. The assumption 
made in previous papers on the origin of the large tran- 
sition noise in YBCO-like and MgB2-like materials, is 
confirmed by the findings of the present work. 



IV. DISCUSSION AND CONCLUSIONS 



The results reported above show several interesting as- 
pects of the transition process in granular superconduc- 
tors with weak and strong links. One aspect, that can 
explain the large noise observed during the resistive tran- 
sition of polycrystalline HTS, is that the transition is 
not a continuous dynamical process. The transition of 
a large number of grains simultaneously occurs to form 
a resistive layer, approximately with the thickness of a 
single grain and orthogonal to the bias current density. 
This permits to evaluate the amplitude of the resistance 
steps generated by the layer formation in real specimens 
on the basis of the average grain size and specimen di- 
mensions. Moreover, a scaling law, deduced from the 
Campbell theorem, permits to deduce the relation be- 
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